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Abstract 

We compute the strange and the average up /down quark masses in the quenched 
approximation of lattice QCD, by using the 0(o)-improved Wilson action and op- 
erators and implementing the non-perturbative renormalization. Our computa- 
tion is performed at four values of the lattice spacing, from which we could ex- 
trapolate to the continuum limit. Our final result for the strange quark mass is 
m^^(2 GeV) = (106 ±2 it 8) MeV. For the average up/down quark mass we obtain 
mf^(2 GeV) = (4.4 ± 0.1 ± 0.4) MeV and the ratio mjm^ = (24.3 ± 0.2 ± 0.6). 



PACS: 14.65Bt, ll.15.Ha, 12.38.Gc. 



1 Introduction 



In recent years, the determination of quark masses has become one of the main research 
topics of lattice QCD simulations. An accurate determination of these masses is in fact 
of great importance for both phenomenological and theoretical studies. The masses of 
the charm and bottom quarks, for instance, enter the theoretical predictions of beauty 
hadron decay rates which, in turn, are relevant for the phenomenological analysis of the 
CKM unitarity triangle and thus of CP violation in the Standard Model. On the more 
theoretical side, a precise knowledge of quark masses may give insight on the physics of 
flavour, by revealing relations between masses and mixing angles or specific textures of 
the quark mass matrices, eventually due to still uncovered flavour symmetries. 

The values of quark masses cannot be directly measured in the experiments, since 
quarks are confined inside the hadrons. On the other hand, being fundamental parameters 
of the theory, quark masses cannot even be computed on the basis of purely theoretical 
considerations. Their evaluation is based on the comparison of the theoretical estimate of 
a physical quantity, which depends on quark masses, and its experimental value. This is 
typically realized on the lattice by using, as experimental input, the values of pseudoscalar 
or vector meson masses. 

In this paper we present the results of an extensive lattice calculation of the strange and 
the average up/down quark masses in the quenched approximation. Particular attention 
has been dedicated to the reduction and control of the systematic uncertainties, particu- 
larly in the case of the strange quark mass. Leading 0{a) discretization effects (where a 
is the lattice spacing) have been removed by using the non-perturbatively (9(a)-improved 
Wilson action and operators [Q]. The systematic uncertainty related to the evaluation 
of the quark mass renormalization constant in perturbation theory has been significantly 
reduced, by implementing the non-perturbative renormalization technique of ref. in 
the RI/MOM scheme. Conversion of the quark masses from the RI/MOM scheme to the 
most popular MS scheme has been performed by using continuum perturbation theory at 
the N^LO 0]. Finite volume effects have also been studied. In order to estimate residual 
systematic uncertainties, we have compared the results obtained by using two alternative 
definitions of the lattice bare quark mass, related to the vector and axial-vector Ward 
identities respectively. Finally, with respect to previous determinations of light quark 
masses by our collaboration [Q, |^, in this study we have performed a calculation at four 
different values of the lattice scale, corresponding to an inverse lattice spacing in the range 
between approximately 2 and 4 GeV. In this way, we have been able to extrapolate our 
results to the continuum limit. The main relevant source of uncertainty, which is left in 
our calculation of the strange quark mass, is therefore the quenched approximation. 

An additional uncertainty is present in the determination of the average up/down 
quark mass. Typical values of the lightest quark masses, used in the present and most 
of current lattice calculations, are approximately of the order of ms/2, where rUg is the 
strange quark mass. Therefore, a large chiral extrapolation is required to reach the phys- 
ical values of the up and down quark masses. Chiral perturbation theory may be used 
as a guidance in this extrapolation, but the inclusion of higher order terms in the chiral 
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expansion, necessary to increase the accuracy of this determination, requires simulations 
with many more (and preferably lighter) quark masses . In the region of masses considered 
in this paper, the pseudoscalar meson mass squared shows a good linear dependence on 
quark masses, and a linear or quadratic fit has been considered in performing the chiral 
extrapolation. The systematic uncertainty introduced by this extrapolation is difficult 
to be reliably quantified. It will be assessed only when simulations on larger lattice vol- 
umes and smaller values of quark masses become feasible. We stress, however, that this 
uncertainty does not affect the determination of the strange quark mass. 

We conclude this section by summarizing the main results of this paper. For the 
strange quark mass and the average value of the up/down quark masses, = (m^ + 
mrf)/2, quoted in the MS scheme at the renormalization scale fi = 2 GeV, we obtain 

mf^(2 GeV) = (106 ± 2 ± 8) MeV (1) 

and 

mf ^(2 GeV) = (4.4 ±0.1 ±0.4) MeV, (2) 

in good agreement with the current lattice world averages 0. For the ratio of the 
strange to the average light quark mass we find 

777- 

— = 24.3 ±0.2 ±0.6. (3) 

This result is in good agreement with the prediction ms/m£ = 24.4 ± 1.5 based on NLO 
chiral perturbation theory [Q. 



2 Details of the lattice calculation 

In determining the values of quark masses we used the standard procedures based on the 
vector and axial- vector Ward identities 0. 

Vector Ward Identity (VWI): in the renormalized continuum theory the VWI reads 

{d,V,ix)O\0)) = (mi(/x) - 7772(^))(5(/i;x)Ot(0)) , (4) 

where 77712 are the quark masses, = qij^q2 is the vector current, and S = qiq2 
is the scalar density. The renormalized quark mass, mg{^) = Zm(/i a) 7^7^(0), is ob- 
tained from the bare mass mq{a) which, on the lattice with Wilson fermions, is equal 
to mq{a) = (l/2a)(l/Kg — l//tcr)- is the Wilson hopping parameter and Hcr is its 
critical value, corresponding to the chiral limit. The VWI relates the quark mass renor- 
malization constant to the one of the scalar density, i.e. = Zg^. The computation of 
the quark mass, using the VWI, can be summarized by the following formula: 

iirS^^\li) = Zm{lia)m!q^^\a) = Zg^{fia) (^1 + bmam^mq{a) . (5) 
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Notice that the parameter bm = —bs/2 provides the ehmination of 0{amq) effects |T]]. 

Axial- Vector Ward Identity (AWI): the renormahzed continuum AWI reads 

{d^A^ix)0^m = 2m(/i)(P(/x;x)Ot(0)) (6) 

where = q'j^'y^q is the axial current, P = q'-y^q is the pseudoscalar density, O is a 
generic operator and we have considered quark fields with degenerate masses. The quark 
mass renormalization constant, in this case, is related by the AWI to the renormalization 
constant of the axial and of the pseudoscalar operators, i.e. Zm = Za/Zp. By choosing 
O = P we have: 

X 

(7) 

where the bare axial current is improved at 0{a) as Aj^^^x) = Afj,{x) + acAd^P{x). The 
coefficient {bA — bp) cancels the terms of 0{amq). For the time derivative, we consider the 
symmetric ((!?(a)-improved) form, i.e. dof{t) = {f{t + a) — f{t — a))/2a. The improvement 
coefficients bm, bA — bp and c^, in eqs. (||) and (|^, are only functions of the bare lattice 
coupling g^. 

Complete information about the lattice calculation performed in this study is provided 
in table |l]. We generated C(IOOO) gauge field configurations in the quenched approxima- 
tion at four values of the gauge coupling constant, corresponding to an inverse lattice 
spacing in the range between approximately 2 and 4 GeV. For each value of the lattice 
spacing, quark propagators have been computed at four light values of the bare quark 
mass, by using the non-perturbatively 0(a)-improved Wilson action 0, p!0[ . 



In view of the final extrapolation of the lattice results to the continuum limit, an 
important requirement in this calculation is a precise determination of the lattice spacing 
which corresponds to the different values of the coupling used in this study. While the 
absolute value of the physical scale is affected by a rather large uncertainty in the quenched 
approximation (of the order of 10%), the ratio between two scales can be determined with 
better accuracy. To this purpose, we use the precise determination based on the study 
of the static quark anti-quark potential 0], which in the range 5.7 < (3 < 6.92 can be 
expressed in the form 

In I J = 1.7331 (/3 - 6) - 0.7849 (p - 6? + 0.4428 (/3 - 6)=^ . (8) 

\a-^{(3 = Q)J 

By using this formula, and varying the inverse lattice spacing at the reference point 
(3 = 6.0 in the conservative range between 1.9 and 2.1 GeV, we obtain the estimates of 
the scale given in table |I[ 
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csw H 
xT 
# conf. 

a-i(GeV) 



6.0 



1.769 



6.2 



1.614 



6.4 



1.526 



6.45 



1.509 



16^ X 52 24^ X 64 32^ x 70 32^ x 70 



500 



200 



150 



2.00(10) 2.75(14) 3.63(18) 



100 



3.87(19) 



K2 



«3 



K4 



0.1335 



0.1338 



0.1340 



0.1342 



0.1339 



0.1344 



0.1349 



0.1352 



0.1347 



0.1349 



0.1351 



0.1353 



0.1349 



0.1351 



0.1352 



0.1353 



0.135175(4) 0.135785(2) 0.135734(2) 0.135680(2) 



Za 11 



(^ = 3 GeV) O 
Z^'ifi = 3 GeV) 



CA m m 



bm P 



{bA - bp) m 



0.804(2) 
0.745(3) 
0.598(3) 
-0.038 
-0.709 
0.171 



0.809(2) 
0.692(3) 
0.575(4) 
-0.038 
-0.691 
0.039 



0.824(2) 
0.668(4) 
0.576(4) 
-0.025 
-0.676 
0.013 



0.825(4) 
0.668(7) 
0.579(8) 
-0.023 
-0.673 
0.010 



amp: fit for t G 
amy'- fit for t G 



[11^25] 
[11 ^23] 
[11 ^ 24] 



[12^31] 
[12 ^ 28] 
[13 ^ 29] 



[17^34] 
[17^28] 
[17^32] 



[17^34] 
[17^30] 
[17^31] 



Table 1: Summary of the lattice details and parameters used in this work. We also give the 
values of the inverse lattice spacing, of the critical hopping parameter and of the renormalization 
constants and improvement coefficients (with corresponding references). In addition we supply 
the reader with the fit intervals that have been used for all the correlation functions considered 
in this work. Note that our time counting is 0, . . . , (T — 1). 
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In the same table, we also give the results for the relevant renormalization con- 
stants, Z^, Zs and Zp, in the chiral limit, which have been determined by using the 
non-perturbative renormalization method of ref. 0], in the RI/MOM scheme. The scale 
dependent constants, Zs and Zp, have been computed at the scale /i = 3 GeV, which lies 
in the allowed range Aqcd < fJ^ < tc /a for all the values of the coupling considered in this 
study. For this reason, our non-perturbative results for the quark masses in the RI/MOM 
scheme will be given at the reference scale = 3 GeV. Details of the non-perturbative 
calculation of the renormalization constants have been presented at the "Lattice 2002" 
conference and will be discussed in a separate publication [plj . 

Concerning the values of the improvement coefficients, ca, bm and (6^ — bp), we use 



the non-perturbative determinations of refs. 10, 12, 13|, whose results are collected in 



table |I[ Notice that at f3 = 6.0 we opted for the value of ca obtained in |]T2| (and recently 



confirmed in ||T5|), whose absolute value is significantly smaller than the one obtained in 
ref. fl^. Had we used the value of ca obtained in ref. 0], we would have found values of 
the AWI quark masses, at /3 = 6.0, smaller by approximately 7%. 

In table |], we show the results for the pseudoscalar and vector meson masses, in 
lattice units, obtained by fitting the corresponding correlation functions at zero spatial 
momentum in the time intervals indicated in table |l[ For each value of the hopping 
parameter, we also present in table ^ the corresponding values of the VWI and AWI 
quark masses, defined in eqs. (|^) and (0), renormalized in the RI/MOM scheme at the 
scale fi = 3 GeV. 

In order to get the physical values of quark masses, we follow the usual procedure 0] 
and fit the pseudoscalar meson masses to the following form 



{ampY = Qi (^am^^^^-* + am2^^^''^ + Q2 (^arrii^^^ + am2^^^''^ , (9) 

and similarly for the AWI quark masses. The subscripts 1 and 2 denote the fiavour of 
the two valence quarks in the meson. Since in this study we only considered mesons 
consisting of two degenerate quarks, we always have mi = 1712 in the fits. For the same 
reason, we did not include a term proportional to (mi — 7712)^ in eq. (§). The physical 
values of the average up/down and of the strange quark masses are then obtained by 
substituting the experimental pion and kaon masses on the l.h.s. of eq. (^ and the values 
of the lattice spacing listed in table |T]. Notice that we do not distinguish the up from the 
down quark mass and, by using eq. (^, we can only determine the average of the two, 
i.e. rrii = {rriu + md)/2. The fit of the pseudoscalar meson masses to eq. (||), at /5 = 6.2, 
and the resulting extrapolation to the physical values, is shown in figure |1], for both the 
VWI and AWI quark masses. As can be seen from the figure, the effect of including a 
quadratic term in the chiral extrapolation of quark masses is rather negligible. This is 
true for all values of the lattice spacing considered in this study. 

The results for the strange and the average up/down quark masses, in the RI/MOM 
scheme, at the renormalization scale /i = 3 GeV, are collected in table |^. Also shown in 
the table are the values of the ratio m^/m^ (which is a both scheme and scale independent 
quantity). Note that this quantity can be determined with much better accuracy than 
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/5 




amp 


amy 




^^(AWI) 


6.0 


0.1335 


0.391(1) 


0.524(4) 


0.0602(3) 


0.0672(4) 




0.1338 


0.356(1) 


0.498(6) 


0.0496(2) 


0.0554(3) 




0.1340 


0.331(1) 


0.480(7) 


0.0425(2) 


0.0475(3) 




0.1342 


0.304(1) 


0.462(9) 


0.0353(2) 


0.0396(2) 


6.2 


0.1339 


0.357(1) 


0.443(3) 


0.0722(3) 


0.0781(4) 




0.1344 


0.303(1) 


0.405(4) 


0.0534(2) 


0.0575(3) 




0.1349 


0.243(1) 


0.370(7) 


0.0343(1) 


0.0370(2) 




0.1352 


0.200(1) 


0.351(11) 


0.0228(1) 


0.0247(1) 


6.4 


0.1347 


0.228(1) 


0.306(2) 


0.0416(3) 


0.0440(3) 




0.1349 


0.204(2) 


0.291(2) 


0.0336(2) 


0.0355(2) 




0.1351 


0.178(2) 


0.277(3) 


0.0256(2) 


0.0271(2) 




0.1353 


0.148(2) 


0.266(4) 


0.0176(1) 


0.0186(1) 


6.45 


0.1349 


0.195(2) 


0.272(4) 


0.0314(3) 


0.0332(4) 




0.1351 


0.167(2) 


0.255(5) 


0.0234(2) 


0.0247(3) 




0.1352 


0.152(2) 


0.247(6) 


0.0194(2) 


0.0205(2) 




0.1353 


0.134(3) 


0.240(8) 


0.0154(1) 


0.0162(2) 



Table 2: Pseudoscalar and vector meson masses together with the corresponding VWI and 
AWI quark masses, in lattice units. Quark masses are renormahzed in the RI/MOM scheme 
at the scale fi = 3 GeV. Information about the fit intervals, the values of the renormalization 
constants and of the improvement coefficients can he found in table ||. 
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P=6-2 



(3=6.2 




0.15 



0.15 



0.00 



0.00 







0.02 



0.04 



0.06 



0.08 







0.02 



0.04 0.06 



0.08 




am (3 GeV) 



Figure 1: Quadratic fits of tlie pseudoscalar meson mass squared as a function of tlie 
renormalized VWI (left) and AWI (right) quark masses, at (3 = 6.2. Empty circles 
represent the lattice data, full diamonds show the physical values of the pion and kaon 
masses. 

the mass itself, since statistical and systematic errors largely cancel in the ratio. 
We emphasize that the results for quark masses presented in table |^ are obtained in a 
completely non-perturbative way. 

3 Conversion to the MS scheme and extrapolation to 
the continuum hmit 

We now convert the RI/MOM quark masses, obtained in the previous section, to the 
popular MS scheme, in which the light quark masses are conventionally expressed at the 
scale fj, = 2 GeV. That allows to compare our results to the results obtained by other 
lattice groups and to the ones obtained by using QCD sum rules. It is only at this stage 
of the calculation that we are forced to introduce perturbation theory. This is because 
the MS scheme, being related to dimensional regularization, is defined in perturbation 
theory only. 

The conversion factor from the RI/MOM to the MS scheme is conveniently calculated 
by introducing the renormalization group invariant mass, m^*^\ defined by dividing out 
from the renormalized quark masses the perturbative scale dependence 



RGI 



ci^i(/i) 



m^S(/i) 

C^(/i) 



(10) 
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f3 


6.0 


6.2 


6.4 


6.45 


myw^ [MeV] 


90(1) 


100(1) 


105(3) 


104(3) 


m^wi [MeV] 


101(1) 


108(1) 


111(3) 


110(4) 


m/wi [MeVl 

t L J 


3.50(6) 


4.02(5) 


4.21(13) 


4.24(23) 


m^wi [MeV] 


3.93(6) 


4.33(6) 


4.45(14) 


4.47(24) 


(m,/m,)vwi 


25.7(1) 


24.90(8) 


25.0(2) 


24.5(6) 




25.6(1) 


24.91(8) 


25.0(2) 


24.5(6) 



Table 3: Values of the strange and the average up/down quark masses in the RI/MOM scheme 
at the scale = 3 GeV, as obtained from the VWI and the AWI methods. We also present the 
values of the scheme and scale independent ratio rris/mi. 



The mass m^*^^ is, by definition, both renormahzation scale and renormahzation scheme 
independent. The beta function of QCD and the quark mass anomalous dimension, enter- 
ing the functions c(/i) in eq. (|10]), are known to 4- loop accuracy, in both the RI/MOM ||] 
and the MS schemes [|16l-||18[. From these papers, we extract 



C^^(/i) 



and 



12 
25 



TT 



12 

25 



TT 



1 + 2.34747 + 12.0599 ( + 84.4076 ( 



1+ 1.01413^+ 1.38921 (^]\ 1.09054 f 



1) 



TT V ^ / \ ^ / 

(12) 

by using ra^;' = 4 as the number of active flavours in the range of scales between 2 and 3 
GeV. Then, by using as{mz) = 0.118, we obtain the conversion factor 

^ mp(2 GeV) _ c^(2 GeV) _ 
^ = m^'iS GeV) " c^'{3 GeV) " ^'^^^ ' ^^^^ 

This result has a N'^LO accuracy in continuum perturbation theory. Therefore, we expect 
the perturbative error in our determination of the MS quark masses to be completely 
neghgible. 

We use eq. ([13|) to convert the results for the RI/MOM quark masses, presented in 
table ^, to the MS scheme. The resulting values of m^^(2 GeV) and m^^(2 GeV), as 
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/3 


6.0 


6.2 


6.4 


6.45 


mywi [MeVl 

s L J 


82(1) 


92(1) 


96(2) 


96(3) 


m^wi [MeV] 


92(1) 


99(1) 


102(2) 


101(3) 


m^wi [MeV] 


3.22(5) 


3.69(5) 


3.87(12) 


3.89(21) 


m|^wi [MeV] 


3.61(6) 


3.98(6) 


4.09(13) 


4.11(22) 



Table 4: Values of the strange and the average up/down quark masses in the MS scheme at 
the scale fj, = 2 GeV. 



obtained from the VWI and AWI respectively, are shown in table |. We do not report in 
the table the values of the ratio m^/m^, presented in table ^, since this quantity is scheme 
independent. 

The last step of our calculation is the extrapolation to the continuum limit. As 
discussed in the previous section, we find it convenient for that purpose to fix the relative 
values of the lattice spacing by using eq. (^), and taking as input the central value a~^(/3 = 
6) = 2.0 GeV. Since our results for quark masses are free of leading 0(a)-effects, the first 
term in the extrapolation is of C(a^). For this reason, we extrapolate our data to the 
continuum limit linearly in a^, and obtain the results 

myw^[ MeV] = (102 ± 2) - (2.01 ± 0.25) , 

mj^\ MeV] = (4.20 ± 0.10) - (0.100 ± 0.013) , (14) 
(m./m^)^^^ = (24.25 ± 0.17) + (140 ± 25) , 

from the VWI, and 

mf^\ MeV] = (106 ± 2) - (1.39 ± 0.28) , 

mf^^[ MeV] = (4.35 ± 0.11) - (0.076 ± 0.014) , (15) 

{ms/mi)^^^ = (24.32 ± 0.17) + (127 ± 25) , 

from the AWI, where the lattice spacing must be expressed in units of fm. The illustration 
of the extrapolation of the MS strange quark mass to the continuum limit is shown in 
fig- i 



9 



VWI 




0.0000 0.0025 0.0050 0.0075 0.0100 
2^2, 

a [fm J 

Figure 2: Extrapolation of the MS strange quark mass, at n = 2 GeV, to the continuum 
hmit. 



A pleasant feature of our results is that the quark masses extrapolated to the contin- 
uum as obtained by using either the AWI or the VWI lead to fully consistent determina- 
tions. The (9(a)-improved quark masses, with the lattice spacings used in our simulations, 
are very close to the continuum limit. Although the central values for the quark masses 
obtained at /5 = 6.45 are slightly smaller than the ones obtained at /5 = 6.4, they are 
completely consistent within the statistical errors. 

4 Systematic uncertainties and final results 

The errors quoted with the continuum determination of quark masses, in eqs. (|14D and 
(|I^), are statistical only. In this section, we discuss the systematic uncertainties and 
present our final results. 

The main sources of systematic errors, present in our calculation, are discussed below. 

o Additive renormalization of the VWI quark mass: when using the VWI method, the 
renormalized quark mass is obtained from the bare one by implementing both a multi- 
plicative and an additive renormalization. The latter is defined by the critical value of 
the Wilson hopping parameter, Kcr, which, in this study, has been determined from the 
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vanishing of the two-point correlation function of the divergence of the axial current. An 
equivalent possibility to fix hi^r is to require the pseudoscalar meson mass squared to van- 
ish in the chiral limit. The values of light quark masses are rather sensitive to the precise 
choice of Kcr, and, when using this alternative determination, we obtain the results 

my^^[MeV] = (110 ±2) - (1.17 ± 0.28) , 

myw^[ MeV] = (4.57 ± 0.11) - (0.052 ± 0.016) , (16) 

{mjmif"^^ = (23.94 ± 0.18) + (20 ± 27) , 

to be compared with those given in eq. ([l^) . By combining the two sets of determinations 
in eqs. (|T^) and (|16]), we get the following estimates of the VWI quark masses 

mj^^ = (106 ± 2 ± 4) MeV , 

m^^^ = (4.38 ±0.10 ±0.18) MeV , (17) 

(m./m^)^^^ = (24.10 ±0.17 ±0.15) , 

where the first error is statistical and the second represents the systematic uncertainty 
due to the spread of the two determinations of Hcr- 

The VWI results in eq. (plT]) are in perfect agreement with those obtained in eq. (plSj) by 
using the AWI method, the differences being smaller than 1%. The systematic uncertainty, 
however, is larger in the VWI case. For this reason, we will quote as our final central 
values of the results obtained from the AWI method, whereas the difference between the 
two methods will be included in the systematic error. 

o Determination of the lattice spacing: our estimate of the lattice spacing has been per- 
formed by using eq. (|]), in which the main source of uncertainty comes from the input 
value of the lattice scale at the reference point (3 = 6. Our choice a~^{(3 = 6) = 2.0(1) GeV 
covers, in a rather conservative way, determinations of the scale based on different phys- 
ical quantities, like f^,, nip, itik*, ro, etc., which are not expected to produce the same 
estimate of the scale in the quenched approximation. In order to evaluate the effect of 
this uncertainty on the determination of the light quark masses, we have repeated the 
analysis by using for a~^{j3 = 6) the values 1.9 and 2.1 GeV respectively, and compared 
the results with those given in eqs. (|I^) and (|l^) obtained by using the central value 
a~^(/3 = 6) = 2.0 GeV. In this way, we find that the quark masses vary by approximately 
5% (the mass increases as the lattice spacing increases). 

o Renormalization constants: a reasonable estimate of the systematic uncertainty involved 
in the non-perturbative RI/MOM calculation of the renormalization constants, Za, Zs 
and Zp, can be obtained by comparing the results with those obtained by using the chiral 
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Ward identity method [[T^, |T^, [TPj . In the latter case, only scale independent quantities can 
be computed. In particular, the values of the renormalization constant obtained from 
the two methods are in perfect agreement within the statistical errors, while a systematic 
difference of the order of 5% is observed in the case of the ratio Zs/Zp. We include this 
difference in the systematic uncertainty. 

o Finite volume effects: the spatial extension of the lattices considered in this study is of 
the order of 1.6-1.7 fm, which is expected to be large enough for finite volume effects to be 
well under control. In order to verify this statement, we have performed an independent 
simulation, sX (3 = 6.0, on the volume 24'^ x 64, which corresponds (in physical units) to 
a spatial extension of 2.3 fm. The results for the VWI and AWI strange quark masses, in 
the MS scheme at the scale fi = 2 GeV, which are obtained from the simulation on the 
larger lattice, are mj^^ = 85(1) MeV and m^^^ = 93(3) MeV respectively. These results 
should be compared with the values mj^^ = 82(1) MeV and mf^^ = 92(1) MeV quoted 
in table ^. From this comparison, we get an estimate of finite volume effects which is of 
the order of 2%, which we also account for in our final systematic error. 

o Continuum extrapolation: the extrapolation of quark masses to the continuum limit has 
been performed by considering only the effect of a linear term in a^. One may wonder, 
however, whether higher order discretization effects are indeed negligible, particularly for 
the results obtained on the coarsest lattice. For this reason, we have also performed the 
continuum extrapolation without including the point at /3 = 6.0 in the fit. In this way, 
we find that the results for the strange and the average up/down quark masses decrease 
by 1% and 2% respectively. 

o Perturbative matching: the conversion factor R of eq. (p!3D, which translates the RI/MOM 
quark masses to the masses in the MS scheme, has been computed by using ^^(m^) = 
0.118 and the number = 4 of active flavours in the range of scales between 2 and 3 GeV. 
By working in the quenched approximation, however, we could have equally computed 
this factor by using up = and as from Aq^^ ~ 0.250 GeV. By proceeding in this way, 
we would have obtained R^^^ = 0.944, instead of the result R^'^^ = 0.918 given in eq. (|T3|). 



The difference between the two determinations, which is of the order of 3%, represents 
an intrinsic ambiguity in the quenched approximation. We include this difference in the 
final evaluation of the systematic uncertainty. 

From the continuum results shown in eq. (0), and by adding in quadrature the sys- 
tematic uncertainties discussed above, we finally obtain our best estimates of light quark 
masses, 

mf^(2 GeV) = (106 ± 2 ± 8) MeV , 

mp(2 GeV) = (4.4 ± 0.1 ± 0.4) MeV (18) 

and 

m,/m£ = 24.3 ±0.2 ±0.6 , (19) 
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which have been also quoted in the abstract of this paper. Notice that most of statistical 
and systematic uncertainties cancel in the ratio ms/m^. 

Our results in eqs. (|1^) and (|T9|) are in agreement with the recent extensive lattice 
QCD calculations of the quark masses performed by using Wilson fermions in the quenched 
approximation |2^-[^. They are also in good agreement with the current lattice world 
averages, presented in the reviews 0, 0] and in the 2002 Review of Particle Physics ^ . 
The main feature of the present study is, in our opinion, the special attention dedicated 
to the reduction and control of the systematic uncertainties within the quenched appro- 
ximation. 
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